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ABSTRACT 

Extended Red Emission (ERE) has been detected in many dusty astrophysical 
objects and this raises the question: Is ERE present only in discrete objects or is it 
an observational feature of all dust, i.e. present in the diffuse interstellar medium? 
In order to answer this question, we determined the blue and red intensities of the 
radiation from the diffuse interstellar medium (ISM) and examined the red intensity 
for the presence of an excess above that expected for scattered light. The diffuse ISM 
blue and red intensities were obtained by subtracting the integrated star and galaxy 
intensities from the blue and red measurements made by the Imaging Photopolarimeter 
(IPP) aboard the Pioneer 10 and 11 spacecraft. The unique characteristic of the 
Pioneer measurements is that they were taken outside the zodiacal dust cloud and, 
therefore, are free from zodiacal light. The color of the diffuse ISM was found to be 
redder than the Pioneer intensities. If the diffuse ISM intensities were entirely due to 
scattering from dust (i.e. Diffuse Galactic Light or DGL), the color of the diffuse ISM 
would be bluer than the Pioneer intensities. Finding a redder color implies the presence 
of an excess red intensity. Using a model for the DGL, we found the blue diffuse 
ISM intensity to be entirely attributable to the DGL. The red DGL was calculated 
using the blue diffuse ISM intensities and the approximately invariant color of the 
DGL calculated with the DGL model. Subtracting the calculated red DGL from the 
red diffuse ISM intensities resulted in the detection of an excess red intensity with an 
average value of ~10 Sio(V)g2V • This represents the likely detection of ERE in the 
diffuse ISM since Ha emission cannot account for the strength of this excess and the 
only other known emission process applicable to the diffuse ISM is ERE. Thus, ERE 
appears to be a general characteristic of dust. The correlation between Nhi and ERE 
intensity is (1.43 ± 0.31) xlO~^^ ergs s~^ sr~^ H atom~^ from which the ERE 
photon conversion efficiency was estimated at 10 it 3%. 



1. Introduction 



Extended Red Emission (ERE) is a broad (AA ~ 800 A) emission band with a peak 
wavelength between 6500 A and 8000 A seen in many dusty astrophysical objects (Figure |I|). For 
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many years, the Red Rectangle was the only object known to possess such an emission feature (see 
Schmidt, Cohen, & Margon [1980] for an excellent Red Rectangle spectrum). This all changed 
with the discovery that many reflection nebulae possessed flux levels in the R and I bands in 
excess of that expected from dust scattered starlight ( Witt, Schild, &: Kraiman 1984 ; Witt & 
Schild 1985| , 1986). The spectroscopic confirmation of the excess flux as ERE (Witt & Schild 1988; 
Witt fc Mahn 19891 ; [Witt et al. 1989| ; [Witt fc Boroson 1990]) proved that ERE was a feature of 
many, but not all, refiection nebulae. The identification of the band as emission was strengthened 
by the imaging polarimetry of NGC 7023 by Watkin, Geldhill, & Scarrott (1991), which showed a 
reduction in the R and I polarization where R and I excess fiux existed. This convincingly proved 
that the excess fiux is due to an emission feature and not changes in the dust scattering properties. 
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Fig. 1. — An example of Extended Red Emission is plotted. This spectrum is for the refiection 
nebula NGC 2327 and has had the underlying scattering continuum subtracted (Witt 1988). 



The detection of ERE in other dusty astrophysical objects quickly followed the confirmation 
of ERE in refiection nebulae. Up to the date of this paper, ERE has been detected in the Red 
Rectangle ( Schmidt et al. 1980 ), refiection nebulae (see above), a dark nebula ( Mattila 1979 ; 
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Chlewicki & Laureijs 1987), Galactic cirrus clouds ([Guhathakurta fc Tyson 1989 ; Guhathakurta| 
fc Cutri 199^ ; ^zomoru fc Guhathakurta 1998] ), planetary nebulae (Furton & Witt 1990, 1992), 
H II regions ( perrin &: Sivan 1992 ; Sivan &: Perrin 1993| ), a nova ( [Scott, Evans, fc Rawlings 1994| ), 



the halo of the galaxy M82 (Perrin, Darbon, & Sivan 1995), and the 30 Doradus nebula in the 
Large Magellanic Cloud (Darbon, Perrin, &: Sivan 1998). 



Clues as to the identity of the material that produces ERE are contained in the above 
observations. The first clue comes from the wide variety of objects that show ERE. The material 
must be able to survive in radically different environments: from cold, quiescent environments 
(dark nebulae and Galactic cirrus clouds) to hot, dynamic environments (reflection nebulae, 
planetary nebulae, and H II regions). Second, the material most likely is carbonaceous since ERE 



has been detected in carbon rich planetary nebulae but not oxygen rich planetary nebulae (Furton 
fc Witt 1992[ ). The robustness of carbonaceous dust material is supported by the essentially 
constant C/H values found for sight lines exhibiting different dust characteristics ( Sofia et alj 



1997), implying that carbon is not easily exchanged between the gas and dust phases of the ISM. 



The third clue comes from the spatial distribution of ERE in reflection nebulae (Witt & Schild 
1988t [Witt fc Malin 1989| ; [Witt et al. T989t [Witt &: Boroson 199q ; [Rogers, Heyer, fc Dewdneyl 



1995 ; Lemaire et al. 1996 ) and planetary nebulae ( Furton &: Witt 1990| ). The ERE is strongest 



in regions where H2 is being dissociated, leading to the conclusion that warm atomic hydrogen 
increases the efficiency of ERE luminescence. So, the material that produces ERE must be robust, 
contain carbon, and produce ERE more efficiently in the presence of atomic hydrogen. 

A prime candidate for this material is hydrogenated amorphous carbon (HAG), which was 
first proposed to explain the ERE in the Red Rectangle ( puley 1985| ). With the discovery of ERE 
in reflection nebulae and other dusty objects, the identification of HAC with ERE has strengthened 
(Puley fc Wilhams 1988| ; [Witt fc Schild 1988| ; Puley fc Wilhams 199q ; [Witt fc Furton 1994| ). The 
identification of HAG with ERE is strongly supported by the laboratory work of Furton & Witt 
(1993). They have shown that the low or non-existent photoluminescence of previously annealed 
HAC or pure amorphous carbon can be greatly enhanced by exposure to atomic hydrogen and/or 
ultraviolet radiation. This corresponds to the strengthening of ERE in H2 dissociation regions 
discussed above. Other ERE producing materials have been proposed, such as filmy quenched 
carbonaceous composite ( ^akata et al. 199^ ) and Ceo ( Webster 1993 ). For a discussion of the 



similarities and differences between these materials see Papoular et al. (1996). 

As ERE has been detected in a large range of dusty objects, the question arises: Is ERE 
present in the diffuse interstellar medium (ISM)? If the answer is yes, this would make ERE a 
general characteristic of dust. A positive answer has been claimed by Duley & Whittet (1990), 



who identified the very broad structure (VBS) seen in extinction curves (e.g. van Breda & Whittet 



1981) as due to ERE. Jenniskens (1994) has pointed out that ERE cannot be the cause of the 
VBS for two reasons. First, spectra showing VBS have had the nearby sky subtracted, effectively 
removing any emission from the diffuse ISM. Second, the VBS strength does not depend on the 
size of the aperture used. Hence, ERE has not been detected in the diffuse ISM. As a result, dust 
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models for the diffuse ISM have not used the abihty to produce ERE as a constraint on possible 



dust grain materials (e.g. 


Kim & Martin 1996 


Dwek et al. 1997 


; Li & Greenberg 1997 


)• 



This investigation is aimed at determining whether ERE is present in the diffuse ISM. 
Detecting ERE in the diffuse ISM is much more difficult than doing the same in a discrete 
object. For a discrete object, ERE detection is done by subtracting a nearby sky spectrum from a 
spectrum of the object (e.g. Witt &: Boroson 19901 ). Subtracting a nearby sky spectrum removes 



contributions to the object spectrum from the Earth's atmosphere (airglow), zodiacal light (dust 
scattered sunlight), and Galactic background light (diffuse ISM, faint stars, and galaxies). This 
results in a spectrum with contributions only from the object being studied, and any ERE is 
directly attributable to that object. As the light from the diffuse ISM is part of the sky spectrum, 
this method will not work for it. A different method is required. 

Two of the strongest (and most difficult to model) sources in a sky spectrum are airglow and 
zodiacal hght ( [Toller 1981] ). Both of these sources can be avoided by simply taking observations 
outside the atmosphere (for airglow) and the zodiacal dust cloud (for zodiacal light). Such 
measurements have already been carried out by the Imaging Photopolarimeters (IPP) carried 
aboard both Pioneer 10 and 11 ( Pellicori et al. 1973| ; [Weinberg et al. 1974| ). The IPP measured 



the intensity of almost the entire sky in the blue (437 nm) and the red (644 nm). By using only 
measurements taken when the Pioneer spacecraft were beyond 3.27 AU, contributions from the 



zodiacal light are avoided ( Banner et al. 1974 ). Therefore, the only known sources contributing to 
the IPP measurements are stars, galaxies, and the diffuse ISM. Using photometric star and galaxy 
catalogs, the contribution to the IPP measurements from stars and galaxies can be removed. The 
resulting blue and red intensities are due only to the diffuse ISM. 

While the IPP measurements have given all-sky maps at two wavelengths and not a spectrum, 
this is sufficient to detect ERE in the diffuse ISM. The presence of ERE in reflection nebulae 
was first detected by observing that these objects had red fluxes in excess of that expected from 
dust scattered starlight ( [Witt, Schild, fc Kraiman 1984[ ). The diffuse ISM is a gigantic reflection 



nebula with the Galaxy's starlight scattered by the Galaxy's dust. Therefore, we can use the same 
criterion, excess red flux, to detect ERE in the diffuse ISM. The scattered light in the diffuse ISM 
is termed Diffuse Galactic Light (DGL). The DGL will have a bluer color than the integrated 
starlight because scattering by dust is more efficient at shorter wavelengths. So, if the diffuse ISM 
color (red/blue ratio) is as red as or redder than the integrated starlight and other sources of 
excess red light can be positively excluded, ERE is present. 

Section ^ describes the Pioneer IPP measurements and the construction of the blue and 
red all-sky maps. The compilation of a star and galaxy photometric catalog, complete to 
approximately 20th magnitude, is detailed in section ^. The detection of ERE in the diffuse ISM 
is contained in section |^. Section ^ presents the properties of the ERE in the diffuse ISM. Finally, 
section |^ discusses the implications of our results and summarizes our conclusions. 
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2. Pioneer Data 



One of the instruments onboard the Pioneer 10 and 11 spacecraft was the Imaging 
Photopolarimeter (IFF). The primary objectives of the IFF were to produce blue and red maps of 
the brightness and polarization of the zodiacal dust cloud from 1 to 5 AU, the background light 
outside the zodiacal dust cloud, and Jupiter ( [Fellicori et al. 1973 ). Of these, we were concerned 
with only the all-sky blue and red surface brightness maps taken outside the zodiacal dust cloud. 

The IFF was a Maksutov-type f/3.4 telescope with an aperture of 2.54 cm and a detector 
consisting of a Wollaston prism, multilayer filters, and two dual-channel Bendix channeltrons 



(Fellicori et al. 1973; Weinberg et al. 1974). Simultaneous measurements were made of the 
orthogonal components of the electric field in both the blue and red. The spectral bandpass 
(half-power) was 3950-4850 A for the blue channel and 5900-6900 A for the red channel ( pellicor 



|et al. 197^ ). See subsection |3.1| for more information on the photometric characteristics of 
the IFF. The IFF instantaneous field of view (FOV) was 2? 29 x 2? 29 for the background hght 
measurements. The IFF was mounted on a movable arm and 64 measurements were taken during 
a single rotation of the Fioneer spacecraft. The angle between the arm and the Fioneer spacecraft 
spin axis (look angle = L) was changed in increments of 1?83 to build up a map of the sky. 
The look angle ranged between 29° and 170° (iFelhcori et al. 19731) . The effective FOV of the 
measurements was 2?29 x (2?29 + 5?625sinL), with a maximum of 2?29 x 7?92 when the look angle 
was 90° and a minimum of 2?29 x 3?27 when the look angle was 170°. At each look angle, a 20 
data roll (rotation) measurement cycle was performed with 8 rolls for the background light, 1 for a 
radioisotope-activated phosphor source (^'^C), 1 for offset and dark current levels, and 10 for data 
readout ( [Felhcori et al. 1973| ; [Weinberg et al. 1974| ; [Toller 198l[ ). 

The raw IFF background sky measurements were processed to produce the Fioneer 10/11 
Background Sky data set available from the National Space Science Data Center (NSSDC). 
The details of the processing can be found elsewhere ( Weinberg et al. 1974 ; [Toller 198l| ; 
Weinberg fc Schuerman 1981 ; Schuerman, Giovane, fc Weinberg 1997| ). A brief description 
of the processing follows. First, the data were calibrated using the inflight measurements of 
the radioisotope-activated phosphor source. Second, the FOV center was computed from the 
spacecraft spin axis direction, the look angle, and the clock angle. Third, the contribution from 
bright stars was subtracted using the stars in the Bright Star Catalogue ( Hoffieit &: Warren 1991[ ) 
and stars with my < 8 ( [Toller, Tanabe, fc Weinberg 1987| ) in the Fhotoelectric Catalog ( Blancc| 
et al. 1968; Ochsenbein 197^ ). Fourth, 37 resolved stars were used to determine the time decay 
of the instrument sensitivity and corrections to the telescope pointing. The final error in the 
positions of the FOVs was on the order of 0?15-0?40. The Fioneer 10 red data have abnormally 
high noise, but as we are also using Fioneer 11 data, this did not adversely affect our results. The 
final Fioneer data are expressed in Sio(V)g2V units, the equivalent number of 10th magnitude (V 



band) solar-type stars per square degree. See subsection [3T[ for details of this unit. 



During the cruise portion of the Fioneer 10 and 11 missions, the IFF mapped the background 
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light a number of times. In order to determine the spatial extent of the zodiacal dust cloud, 
Hanner et al. (1974) examined the brightness of two different regions of the sky as seen by the 
IPP when Pioneer 10 was between 2.41 and 4.82 AU. They found that the brightness of these two 
regions stopped changing after Pioneer 10 passed 3.27 AU, making this distance the outermost 
detectable edge of the zodiacal dust cloud. Therefore, all the measurements taken beyond 3.27 AU 
are free from detectable zodiacal light and useful for this investigation. 

On 5 days while Pioneer 10 was between 3.26 and 5.15 AU and on 6 days while Pioneer 11 
was between 4.06 and 4.66 AU, the IPP mapped the background light in the sky. The resolution of 
a map made on a single day is determined by the FOV of the IPP and its overlap with neighboring 
FOVs. Figure ^a gives an example of the pattern of FOVs using the Pioneer 10 measurements from 
day 68 of 1974. One of the FOVs has been shaded to show the overlap of a FOV with neighboring 
FOVs. The resolution of this map is variable, with each parallelogram being a resolution element. 
In order to actually achieve this theoretical resolution, an algorithm must be used to extract the 
information in the overlapping regions. In fact, the resolution of the IPP measurements can be 
increased significantly by using measurements made on different days. Figure |2|b gives the pattern 
of FOVs for the 11 days used in creating the final high-resolution maps (see below). The FOVs 
from different days do not overlap exactly as the Pioneer 10 and 11 missions were launched on 
different trajectories and the spin axis of each spacecraft changed direction slowly as a function of 
distance from the Sun. 




306 308 310 312 314 316 318 320 306 308 310 312 314 316 318 320 

galactic longitude galactic longitude 



Fig. 2. — The pattern of the FOVs are plotted in these two figures. The pattern of FOVs from 
Pioneer 10 on day 68 of 1974 is shown in (a) with points in the center of the FOVs. One FOV 
is shaded to show the regions which overlap neighboring FOVs. The pattern of FOVs from the 
11 days used in constructing the final maps is shown in (b). Due to the different trajectories of 
Pioneer 10 and 11 and the variable spacecraft spin axis orientation, the FOVs from different days 
do not overlap exactly. 
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2.1. Map Generation Algorithm 



The algorithm used to create the final maps is similar to that by Aumann, Fowler, & Melnyk 
(1990). Their algorithm is called the Maximum Correlation Method (MCM) and it was able to 
improve the resolution of IRAS maps by a factor of ~ 6.5, from ~ 4' to ~ 36". Our algorithm 
was similar to MCM and worked in the following manner. An initial guess at the final image 
(zeroth iteration) was taken as a positive flat image with 0?25 x 0?25 pixels. The next iteration 
was calculated from 



k+l 



1 ^ 



Pi 



(1) 



m=l 



where p'-j is the surface brightness of the kth iteration image at pixel coordinates Cm is the 

correction factor for the mth IPP measurement (the surface brightness in particular FOV) which 
includes pij, and the sum was done over the N IPP measurements which include pij. The value of 
Cm was calculated from 



Cn 



In 



1 

Q 



(2) 



where Im is the mth IPP measurement and the sum was done over the Q pixels which are included 
in the mth IPP FOV. The error, o"^-, in p^- was calculated from 



"a = M 



\ 



N 

E 

m=l 



1 



(3) 



With each iteration, the image gives an improved match to the IPP measurements. The outcome 
of this algorithm is to produce an image which describes all 11 days of the Pioneer measurements. 

The 11 days of IPP background light measurements that were used in constructing the final 
high-resolution maps are tabulated in Table |l|. While the 11 days overall possessed usable data, 
a large number of individual measurements were seen to be of poor quality. There are a number 
of sources for the poor quality data: incorrect subtraction of bright stars, scattered sunlight, and 



corrupt data rolls (Toller 1981). The poor quality data were removed from consideration using 
four criteria. First, data contaminated with scattered sunlight (data taken within 70° of the sun 
for Pioneer 10 and within 45° for Pioneer 11) were removed. Second, all data with negative values 
were removed as these were the result of interuptions in the datastream of between the spacecraft 
and ground station. Third, the data were divided into 5° x 5° boxes and data inside each box 
deviating over 3 standard deviations from the average in either their blue measurements, red 
measurements, or red/blue ratio were removed. Fourth, a small number of points were removed 
by visual inspection. Approximately 25% of the IPP measurements were of poor quality. 

The resulting good data were used as the input for the algorithm described above to produce 
the final high-resolution maps. The algorithm was iterated 10 times, until little change was 
detected. The best iteration map to use depended on the region being investigated. For low 
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Table 1. IPP Usable Days 



Spacecraft 


year 


day 






[years] 


[days] 


[AU] 


Pioneer 10 


1972 


354 


3.26 


Pioneer 10 


1973 


149 


4.22 


Pioneer 10 


1973 


237 


4.64 


Pioneer 10 


1973 


279 


4.81 


Pioneer 11 


1974 


57 


3.50 


Pioneer 10 


1974 


68 


5.15 


Pioneer 11 


1974 


106 


3.81 


Pioneer 11 


1974 


148 


4.06 


Pioneer 11 


1974 


178 


4.22 


Pioneer 11 


1974 


236 


4.51 


Pioneer 11 


1974 


267 


4.66 



'^Sun-spacecraft distance, R, taken 
from NSSDC WWW pages. 
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Fig. 3. — The Aitoff projection (galactic longitude of zero in the center) of the Pioneer blue image of 
the sky is displayed. The resolution of the displayed map is 0?5 by 0?5. The large hole corresponds 
to the Sun's location as seen from Pioneer 10/11. The intensity units are Sio(V)g2V • 
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galactic latitude regions where the amplitude of real structure in the maps is much larger than the 
noise amplitude, the 10th iteration gave the best map. For high galactic latitude regions where the 
real structure amplitude is smaller than the noise amplitude, the 1st iteration gave the best map. 

The 10th iteration map for the blue is displayed in Figure ^ The 10th iteration red map is 
similar to the blue map. Typical uncertainties were 2% for the blue and 3% for the red. Figure ^ 
can be compared directly with the blue background as seen by the Hipparcos star mapper. A map 
of this background is presented in Figure 6 of Wicenec &: van Leeuwen (1995) . The comparison 
is good both in overall strength and morphology. From this comparison, the uniqueness of the 
Pioneer maps was quite apparent as the Pioneer blue map lacks the substantial zodiacal light seen 
in the Hipparcos star mapper blue map. 

As we were only concerned with high latitude regions, we will use the 1st iteration map for 
the rest of this paper and save the higher iterations for later work. The 1st iteration blue and red 
maps are just smoother versions of the 10th iteration maps. 



3. Photometric Star and Galaxy Counts 

In order for this investigation to succeed, the contribution to the Pioneer blue and red 
measurements from stars and galaxies fainter than my = 6.5 was needed. We have tackled this 
problem by constructing a Master Catalog from three separate catalogs, each complete in a subset 
of the range between 6.5 and ~20th magnitude. Ironically, the stars and galaxies with magnitudes 
between 12 and ~20 (Palomar O & E) have the best photometric data available due to the 
existence of the Automated Plate Scanner Catalog of the Palomar Sky Survey I (APS Catalog, 
Pennington et al. 1993| ). In the magnitude range between 9 and 15 (~V band in the north and 



in the south), the Guide Star Catalog (CSC, [Lasker et al. 1990| ; [Russell et al. 199C ; lenknei 



et al. 1990|) provides data in only one band. For the magnitude range between 6.5 and 9.5, there 



exists no good complete photometric catalog. We have used a combination of catalogs (see §|3^ 
to construct a Not So Bright Star Catalog (NSBS Catalog) to give the best currently available 
positions and magnitudes for stars with magnitudes between 6.5 and 9.5. 



3.1. Transformations Between Photometric Systems 



Underlying the construction of the Master Catalog was the transformation between the 
Palomar blue & red (O & E) magnitudes and Johnson B &; R magnitudes to Pioneer blue & 
red (PB & PR) magnitudes. The normalized response curves, R{X), for the blue and red bands 
of all three photometric systems are shown in Figure ^ (Minkowsi & Abell 1963; Lamia 1982; 
Toller 1981|) . The Palomar blue (O) response function was computed for an airmass of 1.5 ( |IIayes| 
& Latham 1975) in order to reproduce the transformation between the Palomar and Johnson 
systems used in calibrating the APS Catalog ( Humphreys et al. 1991 ). For all 6 above bands as 
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Fig. 4. — The response curves are plotted for the Pioneer blue and red channels (PB & PR, Toller 
1981), the Johnson B and R (Lamia 1982) and the Palomar blue and red (O & E, Minkowsi & 
Abell 1963). 
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well as the Johnson V band, the band's equivalent wavelength (Aeq), equivalent bandpass (AAeq), 
zero magnitude flux (F;s^), and the intensity corresponding to Sio(V)g2V and Sio(V)aov units were 
computed and are tabulated in Table ^. The values of Aeq and AAgq are obtained by 

_ JXR{X)dX 

jR{X)dX ^ ' 



and 

JR{X) dX 

'r 



AAeq = (5) 



respectively. The flux corresponding to a magnitude of zero, F;^, in each band was calculated by 
summing the product of the band's response curve and a calibrated spectrum of a Lyrae ( p?ug. 



White, fc Lockwood 1977 ). The calibrated spectrum of a Lyrae was multiplied by 1.028 before 



use to account for the fact that a Lyrae's V magnitude is 0.03 (HofHeit & Warren 1991). The 
intensity corresponding to one Sio(V)g2V unit and one Sio(V)aov unit was computed by summing 
the product of each band's response curve with the spectrum, set to 10th magnitude in the V 



band, of the sun (ILockwood, Tiig, fc White 1992| ) and a Lyrae ( p:iig, White, fc Lockwood 1977|) 



respectively. One Sio(V)x unit is defined as intensity equivalent to one 10th V magnitude star 
of spectral type X per square degree where X is either G2V or AOV. The intensity in mag/o" 
corresponding to Sio(V)g2V units in the B bands, V band, and R bands is 28.5, 27.8, and 27.5 
mag/n", respectively. 

The transformations from the Palomar and Johnson systems to the Pioneer system were 
accomplished by means of the above band response functions (Figure §) and zero magnitude fluxes 
(Table ^) along with an observational grid of stellar spectra spanning the Hertzsprung-Russell 
diagram ( ^ilva &: Cornell 1992 ). This grid consists of spectra covering 3510-8930 A with a 



resolution of 11 A and includes 72 spectral types spanning spectral classes 0-M and luminosity 
classes I-V. Most of the spectra are for solar metallicity stars, but some are for metal-rich and 
metal-poor stars. The spectra were dereddened and stars of similar spectral types were averaged 
to produce the final 72 spectral type spectra ( Silva fc Cornell 1992| ). 



In order to check the accuracy of our transformations, the transformation from the Johnson 
system to the Palomar system was computed and compared to the same transformation as 
determined by Humphreys et al. (1991). Figure |5|a displays the (O—B) correction as a function of 
(B — V) which transforms the Johnson B magnitude to the corresponding Palomar O magnitude. 
Figure ^3 displays the (E — R) correction as a function of (V — R) which transforms the Johnson 
R magnitude to the corresponding Palomar E magnitude. The agreement between the {O — B) 
and (E — R) corrections derived in this paper and those of Humphreys et al. (1991), validates this 
method for deriving transformations between photometric systems. 

The transformation from the Palomar to the Pioneer system is displayed in Figure ^ and 
the transformation from the Johnson to the Pioneer system is shown in Figure |^. We have 
fitted the resulting curves with polynomial functions in order to have an analytic form for the 
transformations. The number of terms in the fitted polynomial was determined by adding terms 
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Tablo 2. Photometric Band Details 



System 


Band 


Acq 


AAeq 


Fa" 




Sio(V)g2V 


Sio(V)aov 






[A] 


[A] 


[ergs cm~^ s~-^ 




[ergs cm" 


"2 s- 


-1 A-i sr-i] 


Johnson 


B 


4467 


1014 


6.632x10" 


-9 


1.198x10" 


-9 


2.174x10-9 


Palomar 





4249 


1168 


6.343x10" 


-9 


1.087x10" 


-9 


2.080x10-^^ 


Pioneer 


PB 


4370 


826 


6.997x10" 


-9 


1.192x10" 


-9 


2.294x10"^^ 


Johnson 


V 


5553 


881 


3.639x10" 


-9 


1.193x10" 


-9 


1.193x10-'^ 


Johnson 


R 


6926 


2057 


1.950x10- 


-9 


8.813x10- 


-10 


6.394x10-1° 


Palomar 


E 


6412 


386 


2.289x10" 


-9 


9.828x10" 


-10 


7.505x10-1° 


Pioneer 


PR 


6441 


968 


2.305x10" 


-9 


9.919x10" 


-10 


7.558x10-1° 



Fa is the flux corresponding to a magnitude of zero. See text for details. 
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Fig. 5. — The transformation from the Johnson system to the Palomar system is plotted. The 
(O—B) correction is displayed in (a) and the {E—R) correction is displayed in (b). Note that the 
iO — B) and {E—R) corrections derived in this paper agree quite well with those from Humphreys 
et al. (1991). 
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Fig. 6. — The transformation from the Palomar system to the Pioneer system is plotted. The 
{PB — O) correction is displayed in (a) and the {PR — E) correction is displayed in (b) . The 
maximum corrections for both [PB — O) and (PR—E) are small. 
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Fig. 7. — The transformation from the Johnson system to the Pioneer system is plotted. The 
(PB — B) correction is displayed in (a) and the (PR—R) correction is displayed in (b). While the 
maximum correction for {PB — B) is small, the maximum correction for [PR — R) is large due to 
the significantly different values of Aeq for the PR and R response curves. 
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until the resulting fitted polynomial followed the general trend of the points. The maximum 
corrections for {PB-O), (PR-E), (PB-B), and {PR-R) are ~0.25, ~0.10, ~0.20, and ~1.0, 
respectively. The large maximum correction for (PR—R) is due to the large difference in the Aeq 
value between the Pioneer (6441 A) and Johnson (6926 A) systems. 



3.2. Master Catalog Construction 

Three star and galaxy catalogs were used in constructing the Master Catalog. The three 
catalogs were the Not So Bright Star Catalog (NSBS Catalog, see below), the GSC, and the APS 
Catalog. As the Master Catalog includes a large number of objects, we chose to construct it in 
small pieces. The GSC is split into 9537 regions ( |Jenkner et al. 1990 ) and we used the same 
regions for our Master Catalog. 

For stars and galaxies with PB magnitudes between 13 and ~20th magnitude, the APS 
Catalog of the Palomar Sky Survey I (POSS I) was used. This ambitious survey is using POSS I 
plates to produce a catalog of stars and galaxies with positions and Palomar O and E magnitudes 
between 12th magnitude and the plate limit (~20th magnitude). Due to difficulties in automated 
photometry in crowded fields, the APS Catalog is limited to |6| > 20° ( [Penning ton et al. 199^ ). 
The APS Catalog is also limited by the lack of any POSS I plates below a declination —33° 



(Minkowsi & Abell 1963). The spatial extent of the APS Catalog defines the regions where we 
can accurately remove the integrated contributions of stars and galaxies to the Pioneer data. The 
positional accuracy is O'.'S ( [Pennington et al. 1993| ). The O and E magnitudes were calibrated using 



UBV photoelectric photometry as detailed in Humphreys (1991). The resulting random error in 
magnitudes was 0.2 for stars and 0.3-0.5 for galaxies ([Pennington et al. 1993| ). The star/galaxy 
classification was done using a neural network ( Odewahn et al. 1992| , 1993). Regions surrounding 



bright stars are not included in the APS Catalog due to scattered light from the bright stars. 
We identified these holes by hand and stars and galaxies from a nearby region of equal size were 
copied into the hole. 

The Guide Star Catalog (GSC) was used to provide data for stars and galaxies with PB 
magnitudes between 9.5 and 13. For each object, the GSC contains a position and a V ((A) w 5600 
A) or a J ((A) = 4500 A) magnitude ( [Lasker et al. 1990| ). The position error is 0'.'2 - O'.'S and 



the magnitude error is ~0.30 mag ( [Russell et al. 1990 ). As the GSC gives only one magnitude. 



deriving accurate PB and PR magnitudes for each object was not possible. Since the GSC is 
complete to 15th magnitude, a large number of the GSC objects also appear in the APS catalog. 
This fact allowed us to statistically transform the GSC object magnitudes to the Pioneer system. 
The average (PB — V) or {PB — J) and (PB — PR) color for stars in common to both the GSC 
and APS Catalog was computed and used to transform the magnitudes of the GSC objects to the 
Pioneer system. 



The catalog of stars between 6.5 V magnitude and 9.5 PB magnitude was constructed 
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by combining three different existing catalogs: tlie Catalogue of Stellar Identifications (CSI, 
Pchsenbein 198^ ), the Michigan Catalogue of Two-Dimensional Spectral Types for HD Stars 
(Houk Catalog, Pouk fc Cowley 1975| ; [Houk 1978| ; [Houk 1982| ; Pouk fc Smith-Moore 1988]) , and 



the UBV Photoelectric Photometry Catalogue (UBV Catalog, |MermilliodT987| |MerminiodT994[ ) 



This produced the Not So Bright Star Catalog (NSBS Catalog) with each star possessing the most 
accurate position, spectral type, PB magnitude, and PR magnitude possible. The CSI was used as 
the base for the NSBS Catalog. It is complete to my ~ 9.5 ( pchsenbein, Bischoff, &: Egret 198l|) 
and is the combination of many other catalogs, notably the Henry Draper Catalog (HD Catalog, 
Cannon fc Pickering 1918-1924| ; pannon 1925-1936| ; pannon fc Walton Mayall 1949|) and the 



Smithsonian Astrophysical Observatory Catalogue (SAO Catalog, SAO Staff 1966| ). Stars already 
subtracted from the Pioneer data were excluded from the NSBS Catalog by excluding all stars in 
the Bright Star Catalogue ( [Hoffieit Sz Warren 199l| ) and stars with mv < 8 ( [Toller, Tanabe,~& 



Weinberg 1987| ) in the Photoelectric Catalog ( [Blanco et al. 19681 ; Pchsenbein 197^ ). 



Information for each star in the NSBS Catalog was gathered in the following manner. The 
star's position was taken from the CSI. The star's spectral type was taken (in order of preference) 
from the Houk Catalog or the CSI. The star's PB and PR photometry was determined from the 
star's Johnson B and R magnitudes which were calculated using one of the following algorithms. If 
the star's V magnitude and (B — V) color were available from the UBV Catalog, the R magnitude 
was computed by using the fit given in Figure ^ to calculate the star's {V — R) color. If only the 
star's V magnitude was available, then the star's (B — V) color was computed from 

{B-V) = iB-V)o + E{B-V) (6) 

where {B — V)o is the star's unreddened color and E{B — V) is the reddening due to dust. A 
similar equation was used to compute the star's (V — R) color. 

The unreddened {B — V)o and {V — R)o were determined using fits between spectral type and 
unreddened color computed using the Johnson response curves and the grid of stellar spectra (Silva 



fe Cornell 1992). Figure ^ displays the fits to the {B — V)o and {V — R)o colors for luminosity class 



V. Fits for luminosity classes HI and I were similar. Ideally, all the stars in the NSBS Catalog 
would have two-dimensional MK spectral types resulting in accurate unreddened colors. This is 



true for stars in the HD Catalog (complete to mpg ~ 9 [ Roman fc Warren 1985 1) with declinations 
less than —12° (Houk Catalog) and a few other stars in the CSI with spectral types from other 
sources ( pchsenbein 198"3| ). For the stars with only one-dimensional temperature classes, the 
unreddened {B — V)o and {V — R)o colors were calculated by taking a weighted average of the 
{B — V)o and {V — R)o colors corresponding to luminosity classes V, HI, and I. The weights were 
the probability of finding a V, HI, or I luminosity class star of the star's temperature class with the 
star's V magnitude and absolute value of galactic latitude. These probabilities were determined 
from the Houk Catalog. 
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Fig. 8. — The transformation between the Johnson (B — V) to (V — R) color is plotted above. The 
open circles were calculated by summing the product of the appropriate response curves and the 
grid of stellar spectra (Silva & Cornell 1992). 
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Fig. 9. — The relationship between a luminosity class V star's temperature class and (a) (B — V) and 
(b) (V — R). The spectral type code (spcode) is a numerical representation of the star's temperature 
class with O = 0, B = 10, A = 20, etc. and subclasses worth their values. For example, a A4 
temperature class has a spcode = 24. 



For each star, E(B — V) and E{V — R) values were calculated from 

-E{B-vy 



and 



E{B-V) 



E{V-R) 



. Av 
E{V-R) 



0.33Ay 



(7) 



A 



V 



0.17yl 



(8) 



respectively ( Whittet 1992|) . The values for E[B — V)/Ay and E{V — R)/Ay correspond to an 
average Milky Way extinction curve, Rv = 3.05 ([Whittet 1992|) . Ay is calculated usmg 



Av = ad 



(9) 



where a = 0.6 mag kpc~^ and d is the star's distance. In order to determine the value of a given 
above, we extracted all 19,746 stars in the NSBS Catalog with two-dimensional spectral types and 
an observed {B — V). Using these stars' two-dimensional spectral types, we calculated their {B — V) 
colors from the algorithm described above with a range of a values. For an a = 0.6 mag kpc"^, 
the calculated and observed {B — V) colors were equal with no galactic latitude dependence of a 
necessary. The distance was computed by solving 



my — My = 5 log 



d 



10 pc 



+ ad 



(10) 



where my is the observed V magnitude and My is the absolute magnitude appropriate for the 
star's spectral type ( Schmidt-Kaler 19821) . Again, for stars without two-dimensional spectral 
types, the computed E{B — V) and E{V — R) values were a weighted average of color excesses 
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Fig. 10. — The combined star/galaxy counts for the NSBS Catalog, the GSC, and the APS Catalog 
are displayed for two Master Catalog regions. A region at b = 19?6 (GSC #2947) is displayed in (a) 
and a region at 6 = 74?6 (GSC #2526) in (b). The counts were normalized to reflect the counts for 
one square degree and one magnitude wide bins. The error bars were computed assuming Poisson 
statistics. 



for V, III, and I luminosity classes. Finally, the star's B and R magnitudes were transformed to 



Pioneer PB and PR magnitudes using the transformations given in §3.1 



We were only interested in the integrated star / galaxy intensity and chose to do the integration 
over 5° x 5° sized regions (see section Since the integration was done over regions larger than 
an individual Master Catalog region which contains ~20,000 objects, the error in the integrated 
intensity was dominated by systematic errors. The random error in the integrated intensity due to 
random errors in the fluxes of individual objects was very small due to the large number of objects 
contributing to the integrated intensity. 

The only systematic error found was associated with the lack of two-dimensional spectral 
types for many of the stars in the NSBS Catalog. By using the population statistics from the Houk 
Catalog (see above), the magnitude of this error was greatly reduced. Using the Master Catalog 
regions with spectral types from the Houk Catalog, we determined the remaining systematic error. 
The integrated flux for the stars in these Master Catalog regions with PB magnitudes between 
5.5 and 9.5 was determined two ways, using the stars' two-dimensional spectral types and using 
only the stars' one-dimensional spectral types. On average, the integrated PB star flux using the 
one-dimensional spectral types was 2.5 it 1.5% too high and the integrated PR flux was 3.3 it 2.5% 
too low. Since there was a systematic zero point error, the PB and PR magnitudes for stars 
without two-dimensional spectral types in the NSBS Catalog were corrected by adding 0.0273 and 
-0.0349 to their PB and PR magnitudes, respectively. This left a random error of 1.5% and 2.5% 
in the integrated PB and PR star flux from the NSBS Catalog, respectively. 
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Fig. 11. — The contribution to the total intensity from each half magnitude bin is plotted for the 
same two regions as in the previous figure. The region at 5 = 19? 6 is displayed in (a) and the region 
at 5 = 74? 6 in (b). The peak contribution comes from about 10th magnitude and the contribution 
from objects with magnitudes > 20 is negligible. 



In order to show that the construction of the Master Catalog worked, we display the 
star/galaxy counts for two Master Catalog regions in Figure 0. Fi ffure PlQa shows the counts for a 
low latitude region and Figure 10 b shows the counts for a high latitude region. These two figures 
show the counts from the NSBS Catalog, the CSC, and the APS Catalog in the regions where 
these catalogs are stated to be valid. The overlap of the three catalogs is good and within the 
expected uncertainties. Additional confirmation of the Master Catalog comes from comparison of 
our star counts with that of the SKY model (Cohen 1994, 1995). The star counts predicted from 
the SKY model agree within Poisson statistics with our observed Master Catalog counts ( Pohen 
1997| ). See Toher (1981) for an excellent review of previous star count work. 



Figure ^ plots the contribution to the total intensity from each half magnitude bin for the 
same regions displayed in Figure |l^. On average, the maximum contribution comes from -^lOth 
magnitude objects. The scatter in the curve for magnitudes < 10 is due to the relatively small 
number of bright stars per Master Catalog region. The contribution from stars and galaxies with 
magnitudes > 20 is negligible and not the result of the Master Catalog incompleteness as the 
curve falls smoothly at magnitudes well below 20. 



4. Detection of ERE in the Diffuse ISM 

The blue and red intensities of the diffuse ISM were determined by subtracting the integrated 
star/galaxy light (ISGL) from the Pioneer maps. As the ISGL was computed from the Master 
Catalog, the spatial extent of the diffuse ISM intensities was limited to regions of the sky covered 
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by the APS Catalog. Currently, the APS Catalog includes approximately 400 plates from the 
POSS I in a fairly random pattern across the sky. Due to the low resolution of the Pioneer 
maps and the size of an individual POSS I plate (6?6, Minkowsi fc Abell 1963| ), regions with 



many adjoining plates were needed. Currently, only two regions have enough contiguous plates to 
accurately derive the diffuse ISM intensities. Region 1 is located between = (355°, 30°) and 
(10°, 60°) and region 2 is located between (/, b) = (90°, -55°) and (120°, -20°). These two regions 
have areas of 315 0° and 820 0° and contain over 110,000 and 746,000 objects, respectively. 

In these two regions, maps of the blue and red intensities for the ISGL were computed using 
the Master Catalog. The fluxes of the stars and galaxies in the Master Catalog were mapped onto 
a grid in galactic longitude and latitude with 0?25 x 0?25 sized pixels. The fluxes were converted to 
Sio(V)g2V units using the conversions given in Table |2| and solid angles (in steradians) computed 
from 

Jl = (/2-/i)(sin62-sin6i) (11) 

where I is galactic longitude (in radians), b is galactic latitude (in radians), and the subscripts 1 
and 2 refer to the minimum and maximum values, respectively. Equation ^ is valid for rectangles 
in galactic longitude and latitude. The resulting ISGL map was smoothed to a resolution of 2° to 
match the Pioneer maps. 

In computing the ISGL intensities, we did not include a weighting in accordance with the 
dwell time of each star within the FOV of individual IPP measurements for the stars fainter than 
m = 6.5. The corrections for the brighter stars (m < 6.5) which were subtracted in the original 
reduction of the IPP measurements (§^, did include such a correction. The reason for this 
difference in treatment was based on the expected surface density of stars of different apparent 
brightness. The fainter stars are sufficiently numerous that many are present in the field at any 
one time, and on average, a star is leaving the field when another similar star is entering the field. 
The passages of brighter stars through the field, on the other hand, are sufficiently singular events 
that the details of these passages must be taken into account. 

Two different cuts in galactic longitude of the Pioneer and ISGL intensities are displayed in 



Figure 12. The first cut was taken from region 1 between galactic longitudes 0° and 5°. The second 



cut was taken from region 2 between galactic longitudes 95° and 100°. Other cuts in regions 1 



and 2 were examined and found to be similar to those displayed in Figure 12. The corresponding 



diffuse ISM intensities were determined by subtracting the ISGL intensities from the Pioneer 



intensities and are plotted in Figure 13. Figure 14 displays the PR/PB ratios (in Sio(V)g2V units) 



for the Pioneer, ISGL, and diffuse ISM intensities. From Figure |T^, it is obvious that the diffuse 
ISM is redder (larger PR/PB ratio) than either the Pioneer measurements or the ISGL. As the 
scattered component of the diffuse ISM (DGL) is bluer (see ^.l| ) than the Pioneer measurements, 
this requires that a nonscattered component is present in the red diffuse ISM intensity. 
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Fig. 12. — The red and blue intensities for the Pioneer measurements and integrated star/galaxy 
hght (ISGL) are plotted as a function of galactic latitude for cuts in galactic longitude between 0° 
and 5° (a) and 95° and 100° (b). Each point corresponds to a 5° x 5° region. The Pioneer error 
bars were computed using the algorithm described in §2.1. The ISGL error bars were assumed to 
be a conservative 5% (see §3.2). 




Fig. 13. — The red and blue intensities for the light from the diffuse ISM are plotted for cuts in 
galactic longitude between 0° and 5° (a) and 95° and 100° (b). The diffuse ISM intensities were 
computed by subtracting the ISGL from the Pioneer measurements. 
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Fig. 14. — Plotted is the red/blue ratio for the Pioneer measurements, the ISGL, and the diffuse 
ISM. The first plot (a) displays the cut in galactic longitude between 0° and 5° and the second plot 
(b) the cut between 95° and 100°. 



4.1. Diffuse Galactic Light Model 



In order to determine the strength of the nonscattered component of the diffuse ISM red 
intensity, the scattered component in the red (DGL) must be removed. An accurate calculation 
of the DGL should include the effects of multiple scattering, the cloudiness of the interstellar 
medium, and the observed anisotropy of the illuminating radiation field. Such a model exists and 



it is the Witt-Petersohn DGL model (WP model, |Witt k Petersohn 1994| ; [Friedmann 1996| ; [Witt 



Friedmann, fc Sasseen 1997| ). The WP model treats the Galaxy as a gigantic reflection nebula. 



The details of the WP model and its use in the ultraviolet can be found in Friedmann (1996) and 
Witt et al. (1997). Below, the salient points of the WP model and its adaptation to the optical 
are described. 

Application of the WP model to the present problem was as follows. We derived the radiation 
field at different heights above and below the galactic plane from the Pioneer all-sky intensity 
maps using the work of Mattila (1980a,b). The all-sky intensity maps were made by adding the 
intensities of the bright stars which were removed during the Pioneer data reduction (Weinberg et 



|al. 1974| ) to the Pioneer all-sky maps after interpolating over the hole caused by the Sun's location 
(see Section |2|) . The total dust optical depth along a particular line-of-sight was computed from 
r = C X Nhi where the Nhi values were from the Bell Laboratories H I Survey ( Stark et al. 1992 ) 
and the Parkes H I survey (pleary, Haslam, fc Heiles 1979| ; |Heiles fc ClearyT979|) . The conversion 
constant, C, was 7.58 x 10"^^ cm^ for PB and 4.71 x 10~^^ cm^ for PR. C was computed using the 
average Galactic extinction curve (|Whittet 1992| ) and {Nhi/E{B-V)) = 4.93 x 10^^ cm'^ 



mag 



(Diplas & Savage 1994). The spectrum of cloud sizes and optical depths, except for scaling to the 
PB and PR bands, was the same as that in Witt et al. (1997). 
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The WP model uses Monte Carlo techniques to describe the radiative transfer through dust. 
With these techniques, photons are followed through a dust distribution and their interaction with 
the dust is parameterized by the dust optical depth (r), albedo, and scattering phase function 
asymmetry (g). The optical depth determines where the photon interacts, the albedo gives the 
probability that the photon is scattered from a dust grain, and the scattering phase function gives 
the angle at which the photon scatters. 

In the blue (PB), reasonable values of the dust albedo and g are 0.50-0.70 and 0.60-0.80, 
respectively ([Fitzgerald, Stephens, fc Witt 1976| ; [Toller 1981| ; [Witt et al. 1982[ ; [Witt, SchildTfc 



Kraiman 1984 ; [Witt, Oliveri, &: Schild 1990| ). Putting limits on the red (PR) values of the dust 
albedo and g was more difficult as the objects best suited to studies of the dust albedo and g (i.e. 
reflection nebulae) are just those objects with appreciable ERE ( [Witt, Schild, Sz Kraiman 1984[ ; 
Witt Sz Schild 19*8^ ). Fortunately, there are a few objects without detectable ERE. One such 



object is the Bok globule investigated by Witt et al. (1990), who found that the dust albedo and 
g values in this Bok globule are smaller by ~30% in the red than in the blue. This albedo and g 
difference between the blue and the red is similar to that predicted by dust grain models (Kim & 
Martin 1996| ; Mathis 1996; Zubko, Krelowski, Wegner 1996[ ; p^i fc Greenberg 1997 ), which predict 
either roughly constant or decreasing albedo and g values when moving from the blue to the red. 
Therefore to be conservative, we have adopted red albedo and g values equal to those in the blue. 
This implies that the DGL intensities we computed in the red are the upper limits. The adopted 
dust grain parameters for the model runs are listed in Table ^. 

A detailed comparison of the diffuse ISM intensities with those predicted for the DGL by the 



WP model is presented in Figures |15| &: 16. The blue diffuse ISM intensities in Figure 15a fall well 
within the WP model predictions, confirming that the blue light from the diffuse ISM is entirely 



attributable to the DGL. This is not the case for Figure 15b where the blue diffuse ISM intensities 
generally have the same galactic latitude dependence as the WP model predictions, but fall 
consistently lower. The differences between the WP model predictions and the actual blue diffuse 



Table 3. WP Model Runs 
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Fig. 15. — The blue intensities for the diffuse ISM as well as the WP model predictions are plotted 
for cuts in galactic longitude between 0° and 5° (a) and 95° and 100° (b). 




Fig. 16. — The red intensities for the diffuse ISM as well as the WP model predictions are plotted 
for cuts in galactic longitude between 0° and 5° (a) and 95° and 100° (b). 
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Fig. 17. — Plotted is the red/blue ratio for the Pioneer measurements, the ISGL, the diffuse ISM, 
and all four of the WP model runs. The first plot (a) displays the cut in galactic longitude between 
0° and 5° and the second plot (b) the cut between 95° and 100°. 



ISM intensities are not surprising as the WP model DGL predictions apply to an average value of 
{Nhi /E{B — V)) , average dust grain parameters (albedo and g), and a radiation field derived from 
the intensities seen at the Earth's position in the Galaxy. The radiation field seen by a dust cloud 
located high above (or below) the Galactic plane could be different from what we calculated using 
the Pioneer maps and the model of Mattila (1980a,b). The dust grain properties along particular 
lines-of-sight are known to vary (e.g. Gordon et al. 199^ and references therein; Calzetti et al. 1995| 



Witt, Friedmann, fc Sasseen 19971) . In addition, the Galactic oxygen abundance gradient found by 
Smartt &: Rolleston (1997) implies a dependence in the dust grain properties and {Nhi / E{B — V)) 
with galactic longitude since oxygen is a primary grain component. In fact, significant variations 
in {Nhi/E{B — V)) have been measured by Diplas & Savage (1994). The disagreement between 



observations and model intensities seen in Figure |16| b makes it clear that a Galaxy-wide study of 
the DGL might reveal important evidence about the dependence of dust properties on galactic 
longitude. 

While the above three points limit the applicability of the WP model DGL intensity 
predictions, the color of the DGL is fairly independent of these points because the DGL color is 
determined almost entirely by wavelength dependence of the dust grain properties (albedo and 
optical depth) between the blue and the red. The values of the dust grain properties we used in 
the DGL model were derived from both observations and dust grain models. The dependence 
of the DGL color on only the wavelength dependence of the dust grain properties is illustrated 



in Figure 17 which plots the PR/PB ratio for the Pioneer, ISGL, diffuse ISM, and DGL model 
intensities. The PR/PB ratio (color) of all four of the WP model runs is approximately unity, 
with a slight gradient in galactic latitude. This figure clearly shows the presence of excess red 
intensity over that expected from scattered red photons as the average observed PR/PB ratio of 



-27- 



the diffuse ISM is approximately two. 

Our results for the presence of an excess red intensity in the diffuse ISM are very similar 
to the findings of Guhathakurta & Tyson (1989) and Guhathakurta & Cutri (1994) who found 
that the {B — R) colors of individual IRAS cirrus clouds were 0.5-2 magnitudes redder than that 
expected for scattered light. The indentification of the red color with ERE in cirrus clouds has 
been confirmed by Szomoru & Guhathakurta (1998) through optical spectroscopy of cloud edges. 
From Figure |l^, the color of the diffuse ISM light is 0.5-4.0 times redder than the DGL. This 
corresponds to a {PB—PR) color 0.4-1.5 magnitudes redder than that expected for the DGL. Our 
results cover a significantly larger region of the the sky (1135 □") than the work of Guhathakurta 
& Tyson (1989) and Guhathakurta & Cutri (1994) which was for individual regions up to 1 



4.2. ERE Intensity 

In order to determine the excess red intensity strength, the contribution from the red DGL 
must be subtracted from the red diffuse ISM intensity. As the WP model is not sufficiently 
complex to reproduce the blue DGL exactly, we could not use the WP model predictions of the red 
DGL. On the other hand, the color of the DGL will not be greatly affected by these deficiencies 
(described at the end of § [4.11 ). Therefore, we used the average DGL color predicted by the WP 
model runs (see Figure p!7| ) along with the observed blue diffuse ISM intensity to calculate the 
expected intensity of the DGL in the red. This results in a predicted red DGL which accurately 
reflects the variations of the dust grain properties, {Nhi / E{B — V)) , and radiation field. 

The red DGL intensities calculated using the above prescription were subtracted from the red 
diffuse ISM intensities to yield a lower limit to the red intensity of the nonscattered component of 
the diffuse ISM. These red intensities along with the Hq intensity expected from H recombination 
in the diffuse ISM ( Reynolds 1984| ) for the two cuts displayed in previous figures are plotted in 



Figure |18[ Clearly, Ha emission cannot explain more than a small fraction of the nonscattered 
component of the red diffuse ISM intensity. This same conclusion was reached by Guhathakurta 
& Tyson (1989) in their study of Galactic cirrus clouds and has recently been confirmed through 



optical spectroscopy of cirrus cloud edges ( Szomoru &z Guhathakurta 199^ ). We identify the 



excess red intensity with ERE, as there are no other known red emission processes in the diffuse 
ISM. The ERE strength we determined is a lower limit on the true ERE strength as the red DGL 
intensity was calculated from the WP model color prediction which was an upper limit on the true 
DGL color. 



5. Properties of ERE in the diflPuse ISM 
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Fig. 18. — The ERE intensities are plotted. The ERE intensities were determined by subtracting 
the predicted red DGL intensities from the red diffuse ISM intensities. The expected intensity from 
diffuse ISM Ha emission (Reynolds 1984) is also plotted using the conversion of 1 rayleigh = 0.26 
Sio(V)g2V • The first plot (a) displays the cut in galactic longitude between 0° and 5° and the 
second plot (b) the cut between 95° and 100°. 



5.1. ERE Correlations 

In order to check our identification of the red excess detected in the previous subsection 
with ERE, we sought correlations of the ERE intensity with known tracers of dust: H I column 
densities ( pleary, Haslam, fc Heiles 1979| ; |Heiles fc Cleary 1979| ; [Stark et al. 1992|) and the COBE 



Diffuse Infrared Background Experiment (DIRBE) 100, 140, and 240 /xm intensities ( [Hauser et al. 
1997t [Toller 19971 ). 



The correlations between H I column densities and ERE intensities for regions 1 and 2 are 



displayed in Figure 19. The points in Figures [19[ a and b were taken from the entire areas of 
regions 1 and 2 and include the representative cuts displayed in Section ^. The best fit line 
for region 1 was y = (2.80 ± 2.05) + (0.0140 ± 0.0036)x. The best fit line for region 2 was 
y = (2.01 lb 3.41) + (0.0159 it 0.0066)x. The best fit lines were determined using the IDL function 
LINFIT. The y-intercepts for both lines are consistent, within the uncertainties, with zero and the 
correlations for both regions are equivalent within the uncertainties. Combining the two regions 
yielded a best fit line of y = (2. 67ibl.71) + (0.0145ib0. 0032)3;. This best fit line gives the relationship 
between the ERE energy emitted per H atom and it is (1.43 ± 0.31)xl0^29 gj.gs s"^ A"^ sr'^ H 
atom~^. 

The correlations between the DIRBE 100, 140, and 240 /xm and ERE intensities are plotted 
in Figure Unlike region 1, there was no clear relationship between the DIRBE and ERE 
intensities for region 2. This was troubling as the correlations between the DIRBE and ERE 
intensities for region 1 was quite good. The origin of the poor correlations for region 2 is related 
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Fig. 20. — This figure displays the correlation between the DIRBE 100, 140, and 240 fim intensities 
and ERE intensity for regions 1 (a) and 2 (b). 
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Fig. 21. — This figure displays the correlation between the DIRBE 100, 140, and 240 ^um intensities 
and H I column densities for regions 1 (a) and 2 (b). 



to the fact that the DIRBE data were collected from a satellite in near-Earth space. While the 
DIRBE data have had the majority of the zodiacal dust emission subtracted, there is a noticeable 



residual left. This is illustrated in Figure 21b which shows the correlations between the DIRBE 
intensities and H I column densities for region 2. The peak in the DIRBE intensities around 
Nhi ~ 550 X 10^^ cm~^ is due to residual zodiacal dust emission as the ecliptic plane runs through 
the middle of region 2. 

The ERE intensity was well correlated with dust tracers for region 1. For region 2, the ERE 
strength was correlated with H I column density but not DIRBE intensities. The poor correlation 
between the ERE and DIRBE intensities was explicable when the presence of residual zodiacal 
dust emission was taken into account. Thus, the identification of the red nonscattered excess 
intensity in the diffuse ISM with ERE was greatly strengthened. Additionally, this is the first time 
in which it has been possible to directly correlate the strength of the ERE emission with known 
dust tracers. 



5.2. ERE Photon Efficiency 

The process which produces ERE is thought to be photoluminescence from a disordered solid 
like HAC. Basically, a UV or blue photon excites an electron from the valence to the conduction 
band and the election relaxes back to the conduction band via a number of (mostly non-radiative) 
transitions which might include one radiative transition. When this radiative transition occurs, an 



ERE photon is emitted (Furton 1994). At most, one ERE photon is emitted per absorbed UV 
or blue photon. The photon efficiency is the percent of photons which are absorbed by the ERE 
producing material which cause the emission of ERE photons. 
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We calculated the lower limit on the photon efficiency by assuming all the photons absorbed 
by dust are absorbed by the ERE producing material. The calculation is fairly straightforward, 
but required a couple of assumptions. The calculation begins with the relationship between the 
ERE intensity and H I column density determined in the previous section. The number of ERE 
photons emitted per H atom was calculated from 

ERE photons ... ( hc\ 

47r^AAe ' ' 



H atoms '^'^ \ A 



cq/ 

_]^4 photons 



5.65 X 10"^* (12) 

s H atom 



where 

ergs 



A = 1.43 X 10"^^ 



s A sr H atom' 

AAeq = 968 A, and Aeq = 6441 A (Table |). The number of ERE photons per unit Ty is then 

ERE photons _ ERE photons Nhi E{B-V)^^^^Av 
TV ~ H atoms E{B-V) Ay ' Ty 

= 9.95x10^ (13) 

S Cm^ Ty 

where E{B-V) /Ay = Ry^ = 3.05"^ (|Whittet 1992|) and Nhi/ E{B-V) = 4.93x10^1 cm^^ mag^^ 



QDiplas k Savage 19941) 



The other half of this calculation was to determine the number of photons absorbed per unit 
Ty. The photons absorbed per unit Ty was calculated from 

abs. photons |-5500 A s (1 ~ e^"^-^ 



/ (1 - ax) ^ ^-IisRpdX 

J 912 A Ty 



Ty J 912 A Ty 

f.5500 A 



/ ^ {l-ax)—IisRFdX (rA«l) 

J912 A Ty 



/912 A Ty 

7.41 X 10« (14) 

S Cm^ Ty 



where ax is the wavelength dependent dust albedo ( pordon et al. 1994 and references therein), tx 



is the wavelength dependent dust optical depth for Ry = 3.05 ( Cardelli, Clayton, &: Mathis 1989[ ) 



and IisRF is the intensity of the interstellar radiation field at the Earth in units of ergs cm ^ s ^ 
A~^ ( Witt &: Johnson 1973| ; Mathis, Mezger, &: Panagia 1983 ). This result is valid in the optically 



thin limit (ta <C 1) which is valid for the diffuse ISM far from the galactic plane. Using the results 
calculated above for the number of photons absorbed and emitted, the ERE photon efficiency in 
the diffuse ISM was 14%. This corresponds to an energy efficiency of 5%. 

This calculation is based on two uncertain assumptions. The first was the strength of the 
interstellar radiation field, Ijsrf ( [Mathis, Mezger, Sz Panagia 1983 ). This Ijsrf is likely to be too 



low as the Mathis et al. (1983) Iisrf was calculated for a smooth diffuse ISM, whereas the diffuse 
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Fig. 22. — The energy emitted per H atom by the diffuse ISM is plotted. The ERE point is from 
this work. The COBE points are from the work of Arendt et al. (1997) and Boulanger et al. (1996). 



ISM is known to be clumpy ( Elmegreen fc Falgarone 1996| ) and the radiative transfer in a clumpy 



medium is greatly different from that in a smooth medium ( |Witt &: Gordon 1996| ). Additional 
discussion of the appropriate Ijsrf can be found in Mathis (1997). 

We determined the actual level of the Iisrf by equating the energy absorbed by dust to 
that emitted by dust. Except for ERE, the energy emitted by diffuse ISM dust is emitted in 
the infrared and the Cosmic Background Explorer (COBE, [Boggess et al. 1992 ) contained the 



Far-Infrared Absolute Spectrophotometer (FIRAS) and DIRBE which measured the infrared 
emission from dust. Using data from high latitude regions, the diffuse ISM spectrum per H atom 
from 3.3 to 240 /.im was derived by Arendt et al. (1997, as presented in Dwek et al. [1997]) and 



from 100 yumto 1 mm by Boulanger et al. (1996). In Figure 22, we have plotted the spectrum of 
the diffuse ISM including the contribution from ERE as determined in this work. Integrating this 
spectrum gives the emitted energy which is 5.77x10"^^ ergs s~^ H atom~^. The ERE contributes 
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a quite significant 3% of tlie total energy emitted by dust. In order to get the same amount of 
energy absorbed by dust between 912 A and 1.2 fim, the Mathis et al. (1993) Iisrf must be 
multiplied by 1.33. Using this elevated Iisrf and given the uncertainties of all the quantities 
contributing to this result, we calculated the lower limit ERE photon efficiency to be 10 it 3% and 
the corresponding energy efficiency to be 4 it 1%. 

The second assumption was that all the photons absorbed by dust were absorbed by the 
material which produces ERE. This is unlikely to be true. Reducing the number of photons 
absorbed by the ERE producing material would raise the ERE efficiency. Assuming a 100% 
conversion efficiency for the ERE producing material would require that 10% of the photons 
absorbed by the dust are converted to ERE photons. This makes the ERE producing material an 
important component of dust grains. 

Instead of a wavelength independent ERE conversion efficiency, perhaps it is a question of a 
high efficiency in a narrow range of wavelengths. Such behavior could be caused by a strong dust 
absorption feature such as the 2175 A extinction bump. The percentage of the photons absorbed 
by the 2175 A bump material as compared the total number of photons absorbed by dust, was 
calculated as the difference between the number of photons absorbed for dust between 912 A and 
5500 A with and without a 2175 A bump . The number of photons absorbed by bump less dust was 
calculated using equation 0, a Ry = 3.05 tx curve from Cardelli et al. (1989) with the bump term 
removed, and a\ = 0.6. We calculated that 13% of the photons absorbed by dust are absorbed by 
the 2175 A absorption feature. This percentage is higher than the 10% ERE efficiency and, thus, 
there are enough photons absorbed by the 2175 A bump material to produce the observed ERE 
intensity. We are not claiming the 2175 A bump is the source of the ERE photons, only that this 
is the kind of absorption feature needed if the ERE photons are absorbed in a narrow range of 
wavelengths. It is intriguing that a known absorption feature, which has never been definitively 
identified with a particular material, could be the source of the ERE photons. 



6. Discussion and Conclusions 



The detection of ERE in the diffuse ISM has direct consequences for all dust grain models. 
With this work, ERE has been detected from H II regions to the diffuse ISM effectively spanning 
the entire range of environments where dust is present. Thus, ERE is a characteristic of dust in 
general. None of the current dust grain models (e.g. [Kim &: Martin 1996| ; [Mathis 1996 ; Zubko 



Krelowski, Wegner 1996 ; Dwek et al. 1997 ; Li fc Greenberg 1997 ) specifically includes a material 



which can produce ERE. As ERE has been detected only in carbon rich planetary nebulae ( [Furton 
fc Witt 1992| ), this adds ERE to the already long list of dust features attributed to carbonaceous 
materials. This further deepens the current "carbon crisis" ( ^now fc Witt 1995[ , 1996) - the 
conflict between the small amount of carbon available for dust grains and the large number of 
observed dust features attributed to carbonaceous materials. As the identification of other dust 
features with carbonaceous materials is not as strong as that of ERE, the connection of some dust 



features with carbonaceous materials needs to be reevaluated. 



6.1. Anomalous Dust Scattering Properties instead of ERE? 

A possible (but highly unlikely) explanation of the red excess observed in the diffuse ISM is 
anomalous dust scattering properties in the red. If the dust albedo increased to unity (perfectly 
scattering dust grains) around 6500 A, it would almost be possible to explain the strength of the 
red excess, because the red excess is approximately equal to the scattered intensity which was 
computed for albedos between 0.5 and 0.7. In fact, this argument was presented by Greenstein 
& Oke (1977) in their attempt to explain the excess red intensity of the Red Rectangle. They 
concluded that only with most favorable geometry and a high albedo could they explain the excess. 

A scattering explanation has been effectively ruled out by means of spectropolarimetric 
results, which show that the polarization decreases where ERE increases. For example, this 
explanation for the Red Rectangle was convincingly shown to be false by the spectropolarimetry of 
Schmidt et al. (1980). Schmidt et al. (1980) observed a reduction in the polarization across the red 
excess feature which proved that the excess was due to emission and not scattering. The emission 
nature of the red excess in reflection nebulae has also been confirmed by the imaging polarimetry 
of NGC 7023 by Watkin et al. (1991), who found that the polarization decreased in locations were 
the ERE was the strongest. The two polarimetric observations discussed above are the only two 
polarimetric observations which have been done on objects with ERE. This suggests that the red 
albedo is approximately equal to the blue albedo as derived by Witt et al. (1990) for a Bok globule 
without ERE and calculated from dust grain models ( |Kim fc Martin 1996| ; [Mathis 1996 ; Zubko 



Krelowski, Wegner 1996| ; Li &: Greenberg 1997|) . Thus, anomalous red dust scattering properties 



are extremely unlikely to account for the red excess in the diffuse ISM. 



6.2. Iere/ IscA and ERE efficiency 



The photon efficiency calculated above for the ERE in the diffuse ISM is 10%. The question 
arises: Is this efficiency characteristic of ERE in general or only ERE in the diffuse ISM? It should 
be possible to calculate the ERE efficiency knowing the scattered intensity {Isca), the ERE 
intensity, and the type of illuminating radiation field. From the scattered intensity and the type of 
radiation field, the number of photons absorbed by dust in the UV and blue should be calculable. 
In practice, the relationship between the scattered and absorbed intensity is quite complicated 
due to the non-isotropic nature of scattering by dust grains (see §4.1). 



Figure plots the ERE intensity versus the scattered intensity for all objects with ERE and 
scattered intensities available in the literature ( Mattila 197S| ; Witt &: Boroson 1990; Perrin & 



Si van 199^ ). This figure is similar to Figure 1 in Witt &: Boroson (1990) but with a much larger 



range in Iere and IscA values due to the inclusion of Orion H II region points, at the high end. 
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Fig. 23. — This plot shows the ERE versus scattered intensities for various objects were the ERE 
intensity has been measured. The diffuse ISM points are from this work. The L1780 dark cloud 
point is from Mattila (1979). The reflection nebulae and the Red Rectangle points are from Witt 
& Boroson (1990). The Orion H II region points were taken from Figures 3 &; 4 of Perrin h Sivan 
(1992). The dashed lines show Iere/Isca ratios of 1%, 5%, 20%, 100%, and 200% and are labeled 
appropriately. 
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and diffuse ISM points, at the low end. The Iere/Isca ratio ranges from 0.01 to over 2.0. The 
Red Rectangle, which was long thought to be unique in its ERE intensity, is seen to have normal 
Iere/Isca ratios. 

Interpreting the scatter in the Iere/Isca ratio involves detailed knowledge of the geometry 
and type of illuminating radiation field in each object. The dust scattering geometry affects 
Iere/Isca such that the larger the angle of scattering by dust grains, the lower IscA for a given 
radiation field, and the larger the Iere/Isca value. Essentially, the scattered intensity is reduced 
due to the forward scattering nature of dust grains, while the ERE intensity is unaffected as the 
dust emits ERE isotropically. The illuminating radiation field affects Iere/Isca since the Iere 
intensity is proportional to the radiation field between 912 and 5500 A and IscA is proportional 
to the number of photons between 5500 and 8000 A. Thus, the bluer the illuminating radiation 
field, the larger the Iere/Isca- 

The lack of low values of Iere/Isca in the diffuse ISM (0.05-2.0) is likely due to the fact 
that the diffuse ISM observations were all at high galactic latitudes where the dust predominately 
scatters photons from the Galactic disk at large angles. Similarly, the high value of Iere/Isca 
(0.3) in the dark nebula L1780 is due to its location high above the Galactic disk (b ~ 36°) 
and resulting large scattering angle. The bipolar geometry of the Red Rectangle, oriented 
perpendicular to our line-of-sight, leads to predominately large angle scattering, again giving large 
Iere/Isca values (0.5-1.5). The range in Iere/Isca (0.01-0.68) in reflection nebulae is mostly 
due to the geometry and this is superbly illustrated as the highest reflection nebula Iere/Isca 
value (0.68) is for IC 63, an externally illuminated nebula which scatters photons at ~ 90° (|Witt 



et al. 1989; pordon et al. 19*97 ). The real range of Iere/Isca in reflection nebulae is larger than 



that represented by the spectroscopic detections of ERE shown in Figure 23. Figure |2^ plots the 
(B — V) and (V — R) colors observed in 14 reflection nebulae by Witt & Schild (1986). As can be 
seen from the location of the dashed line, the maximum Iere/Isca is at least 1.0. On the other 
hand, the Orion H II region has high values of Iere/Isca (0.2-0.6) due to the extreme blueness 
of the illuminating radiation fleld which originates from the hot Trapezium stars (late O, early 
B spectral types). Thus, the scatter in Iere/Isca could easily be explained by geometry and 
radiation field effects with a constant ERE efficiency of 10%. The test of this hypothesis awaits 
a detailed calculation of the ERE efficiency in each object properly accounting for the effects of 
geometry and illuminating radiation field. 



6.3. Conclusions 

Using blue and red all-sky measurements taken by Pioneer 10 and 11 outside the zodiacal 
dust cloud along with star and galaxy counts to 20th magnitude, we determined the blue and red 
intensity of the diffuse ISM in two large regions with areas of 315 0° and 820 0°. By comparison 
with a model for the DGL, the blue diffuse ISM intensity was found to be entirely attributable to 
DGL. The diffuse ISM red intensity was found to consist of DGL and ERE in roughly equal parts. 
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Fig. 24. — The A{B — V) and A{B — R) colors of reflection nebulae observed by Witt & Schild 
(1986) are plotted. The A{B — V) color is the nebular (B — V) color minus the illuminating star's 
(B — V). The A{B — R) color is defined similarly. The solid line shows the expected relationship 
between A{B — V) and A{B — R) for scattering assuming the dust albedo is equal at B, V, and R 
wavelengths. The dashed line gives the relationship for Iere/Isca = 1- 
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Thus, the ERE is detected in the diffuse ISM and shown to be a general characteristic of dust in 
all dusty environments. The ERE in the diffuse ISM is consistent with photoluminescence with an 
photon efficiency of 10 it 3% (lower limit) which corresponds to a 4 ± 1% energy efficiency. 

We plan to expand on this work by expanding the fraction of the sky studied. The APS 
Catalog is scheduled for completion by the end of 1997 ( Humphreys 199^ ) and will cover ~50% of 



the sky. We anticipate using the SKY model ( [Cohen 1994 , 1995) to derive the faint star counts in 



regions not covered by the APS catalog after calibrating the SKY model in those regions where 
the APS catalog exists. This will allow us to investigate the distribution of ERE and DGL over 
almost the entire sky. The all-sky distribution of ERE will allow us to determine the global 
ERE photon efficiency as well as any variations as a function of position in the sky. The all-sky 
distribution of the DGL will allow us to quantify the variations in the WP model inputs (dust 
grain albedo and g, {Nhi/E{B — V)), and radiation field) on a Galaxy-wide scale. In addition, the 
comparison between the SKY model and the APS Catalog will yield valuable information about 
Galactic structure. 
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